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ABSTRACT 

In the recent literature there is circumstantial evidence that the viscosity of the 
intracluster medium may not be too far from the Spitzer value. In this letter, we 
present two-dimensional hydrodynamical simulations of ram pressure stripping of disc 
galaxies in a viscous intracluster medium (ICM). The values of viscosity explored 
range between 0.1 and 1.0 times the Spitzer value. We find that viscosity affects the 
appearance and the dimensions of the galactic wakes but has very little effect on the 
evolution of the gas mass of the galaxy. 

Key words: galaxies: spiral - galaxies: evolution - galaxies: ISM - galaxies - indi- 
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1 INTRODUCTION 

In clusters, galaxies can lose some or all of their gas 
by ram pressure stripping (RPS) due to their motion 
through the i ntracluster medium (ICM). Both, analyti- 
cal estimates dGunn fc Gotd 1971) and (hydro)dynamical 
simulations (e.g . lAbadi et all Il999l. iQuilrs fc Moorei [200ll 
ISchulz fc Struck! l200ll, IVollmer et"atl l200ll. iMarcolini et all 
12001 IRoediger fc BruggenllioO^ . IRoediger fc BruggerJ[2007l ) 
show that RPS can remove a significant amount of gas from 
galaxies, and is thus important for the evolution of galaxies 
and the ICM. 

In addition to gas loss by ram pressure pushing, galaxies 
also suffer gas loss by continuous (sometimes also called tur- 
bulen t or viscous) strip p ing (e . g. lNulsenlll982l: [Quilis et all 
2000: ISchulz fc Struckl l200ll: iRoediger fc Hensler 



2005; 



Roedigcr & Briiggcn 2006). 

iNulsenl (|l982l ) has studied the effects of transport pro- 
cesses and turbulence on the flow of gas past a galaxy and 
has found that they could produce more stripping of gas 
than ram pressure alone. For turbulent stripping, he found 
a mass stripping rate of 
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where r = lOrio kpc, picm = 10~ ri-zmn g cm and 
v = 1000«3 km s~ . This is similar to wha t is fo und in the 
simulations by IRoediger fc Briiggenl l|2006l . 120071 ) who find 
rates of about ~ IMqvt -1 . Turbulent stripping is mainly 
caused by the Kelvin-Helmholtz instability, which is sup- 
pressed by viscosity th at stabilises modes with wav elengths 
smaller than r/Re (e.g. iBetchov fc Criminalelll967l ). 
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There exist only few detailed studies on the wakes 
of stripped galaxies. Using the hydrodynamical adaptive 
mesh refinement code FLASH, we have s tudied galactic 
ram pressure tails in a constant ICM wind l)Roediger et al.l 
2006) as well as for gala xies on realistic cluster orbits 
|Roediger fc Briiggenl 120081 ). The minimum tail width of 
about 20 to 30 kpc is found near the galaxy. With in- 
creasing distance to the galaxy, the tail flares to widths 
of 30 to 80 kpc at a distance of ~ 100 kpc behind the 
galaxy. Ot her hydrodynamical simulations using either g rid 



codes (e.g. iQuilis fc Moore] l200ll . iMarco ini etaD l2003h or 



Abadi et alj 1 19991 . 



smoothed particle hy drodynamics (e.g. 
ISchulz fc Struck 2001) do not focus on the tails. However, 
according to the snapshots provided in these papers, the 
gas tails are similar to the ones in our simulations. Using a 
sticky-particle code, Vollmer et al. (e.g. 1999, 2000, 2001a, 
2001b, 2003, 2004a, 2005, 2006) aim at reproducing the RPS 
history of individual galaxies by comparing simulations and 
observations. However, these simulations concentrate on the 
gas distribution close to the galaxy. 



Recently, lOosterloo fc van Gorkoml (|2005l ) presented 
observations of a ~ 120 kpc long and < 25 kpc wide tail 
of H I gas associated with NGC 4388, and also suggest that 
this tail is due to ram-pressure stripping of this galaxy, ei- 
ther in the ICM of the Virgo cluster or in the halo of the 
nearby elliptical galaxy M86. However, the tail of NGC 4388 
is narrower than the tails i n the simulations presented in 
IRoediger fc Briiggenl (|2006l . 120081 ). although it seems to be 
flaring in a similar fashion as found i n the simulation . 
Also the X-ray a nd Ha tails observed bv lSun et al] (2006); 
ISun et all ([20071 ) : I Yagi et al l l|2007h are much narrower than 
simulated ram pressure tails: While being only ~ 7 kpc wide, 
they reach lengths of ~ 70 kpc and show hardly any flaring. 



© 2008 RAS 



2 E. Roediger and M. Briiggen 



This difference may be caused by the microphysics of the 
ICM. 

The viscosity of the intracluster medium has been dis- 
cussed before. Based on observations of the P erseus cluster 
it has been suggested bv lFabian et al.1 (|2003l ). iFabian et al.l 
that viscosity may play an important role in dissi- 
pating energy injected by the central AGN. Circumstan- 
tial evidence for the presence of significant ICM viscosity 
is also provided by an examination of the morphology of 
Ha filaments in the Perseus cluster. Several of the filaments 
appear to trace well-defined arcs which argues against the 
presence of strong turbulence in the ICM core, possibly re- 
sulting from the action of vi scosity. This idea has been tested 
in numerical simulatio ns bv lRuszkowski et al.l (|2004bl lah and 
iRevnolds et"ai1 i|2005l ). 

In the case of a fully ionized and unmagnetised, ther- 
mal plasma, therelevant coefficie n t of v iscosity is given 
by iBraginskiil (|l95Sl ) and ISpitzeij (|l962h as p ps 6.0 x 
10" 17 (lnA/37)" 1 T 5/2 g cm" 1 s" 1 , where T is the temper- 
ature of the plasma measured in Kelvin and In A is the 
Coulomb logarithm. It results from the cumulative effect of 
weak Coulomb collisions. The mean free path of such inter- 
actions scales with T 2 . As the sound speed scales with T 1 ^ 2 , 
the viscosity is proportional to T 5 / 2 . 

It is customary to measure the importance of viscosity 
through the Reynolds number, Re = vl/v, where v and I 
are characteristic velocities and length-scales of the system 
and v = pip the kinematic viscosity, with p being the fluid 
density. For our case 



Re ~ 26t! 3 rion_3/„ 
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which indicates that viscosity may play a role unless the 
viscosity is strongly suppressed, i.e. the suppression factor, 
/„, is sufficiently small. 

This suppression results from the cluster magnetic 
fields. Even weak fields lead to a tiny proton gyroradius, 
which results in a very effici ent suppressi on (factor of ~ 10 23 
for typical ICM conditions, ISrjitzerll 1962h of the local viscos- 
ity perpendicular to the magnetic field. Mag netic fields in 
the ICM are certa inly tangled or even chaotic l|Clarkell200J : 
Enfili rTet alj|2005h and will lead to a reduced macroscopic 
viscosity. The degree of reduction, however, is unknown. Due 
to the same mechanism, also thermal conduction in the ICM 
is suppressed perpendicular to the magneti c field lines. In a 
recent study, iNaravan fc Medvedevl (|200ll ^ found an effec- 
tive macroscopic thermal conductivity that is a factor of 
f v ~ 10 -2 — 0.2 lower than the unmagnetised value. Similar 
arguments may apply to the vis cosity. From stud i es of ICM 
turbulence in the Coma cluster, ISchuecker et al.l l|2004h de- 
rive an upper limit on the kinematic viscosity of the ICM 
of ~ 3 • I0 29 cm 2 s~ 1 . For typical ICM densities and tem- 
peratures, this corresponds to a viscosity suppression factor 
f v around 0.1. Thus, we might expect the viscosity to be 
suppressed by some factor between 10 -2 and unity. 

The magneti c fields of the ICM (e.g. IClarkd 120041 : 
EnBlin ct al. 2005|) may themselves influence the appearance 
of ram pressure stripped galaxies. Given that the thermal 
pressure in the ICM dominates over the magnetic pressure, 
the magnetic fields should be frozen-in and follow the ICM 
flow. Thus, they should be generally parallel to the galaxy's 



tail. Such a magnetic field structure could attenuate ther- 
mal conduction between the stripped gas and the hot ICM, 
and it could suppress e.g. Kelvin-Helmholtz and maybe even 
Rayley- Taylor instabilities in the tail. However, it is unclear 
how well the magnetic fields will be aligned with the ICM- 
ISM interfaces and how strong this suppression will be. The 
influence of the magnetic fields will be studied in a forth- 
coming paper. 

Here, using adaptive-mesh, hydrodynamical simula- 
tions, we investigate the effect of a macroscopic viscosity 
on the stripping of gas from a galaxy. 



2 METHOD 

We model the ICM-ISM interaction in the galaxy's rest 
fr ame, i.e. the galaxy is expos ed to an ICM flow. The work 
of iRoediger fc Briiggen! (|2007l ) showed that the classical es- 
timate of the stripping radius given bv lGunn fc Gottl (|l972T ) 
can be adapted to galaxies exposed to a variable ICM wind. 
Therefore, here we decouple the effect of time- variability and 
use a constant ICM wind to focus on the effect of viscosity. 

To include viscosity in our hydro-simulations, we add 
velocity diffusion to the momentum equation 
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where P is pressure, p density, Vi the components of velocity, 
gi the components of the gravitational acceleration, 5ik the 
Kronecker delta, 



eik = 



1 / duj du k \ 



2 \dx k dxi J 



(•>) 



and A = en. Here p is the coefficient of (shear) viscos- 
ity. In our simulations we neglected bulk viscosity since it 
vanishes for an ideal gas. We use the standard Spitzer vis- 
cosity for an unmagnetised plasma, for which p = 6.0 x 
10~ 17 {lnA/37)' 1 T 5/2 f v g cm" 1 s" 1 . 

We have noted earlier that the precise value of the sup- 
pression factor is highly uncertain and, depending on the 
nature of magnet i c turb ulence, may even exceed the Spitzer 
value (|Cho et al] l|2003h ) or be suppressed well below it. For 
our viscous simulations we choose f v between and 1 (see 
Tabled. 

In order to ensure numerical stability and to prevent 
too small timesteps we have switched off viscosity for tem- 
peratures outside the range 3 ■ 10 5 K < T < 1.5 • 10 s K. 

The total energy fluxes are not modified by viscosity as 
the effects of viscous heating are small. 



2.1 Code 

The simulations w ere performed with the FLASH code 
|Frvxell et aU fcoOO) version 2.5, a multidimensional adap- 
tive mesh refinement hydrodynamics code. It solves the Rie- 
mann problem on a Cartesian grid using the Piecewise- 
Parabolic Method (PPM). 

The viscous runs are fairly expensive because the 
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Table 1. ICM conditions for different simulation runs. 



run: 


M08 


M08-V1 


M08-V2 


M08-V3 


PlCM/10 _27 gcm- 3 


f 


1 


1 


1 


«ICM/1000 kms" 1 


0.8 


0.8 


0.8 


0.8 


fv 





0.1 


0.5 


1.0 



run: 


M20 


M20-V 


PlCM/10 _27 gcm- 3 


1 


1 


^iCm/1000 kms" 1 


2 


2 


fv 





0.1 



timestep imposed by the viscosity scales as (Ax) 2 /fj,, where 
Aa; is the resolution of the computational grid. As this 
timestep scales more strongly with Aa; than the standard 
hydrodynamical Courant condition, and because viscosity 
depends strongly on temperature, the constraints on the 
timestep are more stringent than in the inviscid runs. Hence, 
here we present only results from 2D simulations in cylin- 
drical coordinates, (R,Z). 

The ICM wind enters the simulation box at Z = 
— 65kpc. The boundary at R — is the galaxy's symme- 
try axis, the remaining two boundaries obey outflow condi- 
tions. We use a grid size of 260 kpc in ^-direction 65kpc in 
/indirection. With 5 refinement levels we reach an effective 
number of 2048 x 256 grid cells and effective resolution of 
0.25 kpc. 

2.2 Model galaxy 

The g alaxy model is the same as in iRoedieer fc Briiggenl 
(2006), i.e. a massive spiral with a flat rotation curve at 
200 kms -1 . It consists of a dark matter halo (1.1 ■ 10 11 Mq 
within 23kpc), a stellar bulge (10 10 Mq), a stellar disc 
(10 11 Mq) and a gaseous disc (5 • 10 9 Mq). All non-gaseous 
components just provide the galaxy's potential and are not 
evolved during the simulation. For a description of the in- 
dividual components and a list of parameters please refer 
to RB06. Initially, the gas disc is set in hydrostatic equilib- 
rium with the surrounding ICM (see also RB06). The disc's 
rotation is included via the centrifugal force. 

2.3 ICM conditions 

Table[T]lists the ICM conditions for the simulation runs. The 
ICM temperature was 7.2 • 10 7 K or 6.2 keV in all runs. We 
start the simulation with the ICM at rest and then increase 
the inflow velocity over the first 50Myr from zero to wicm 
(for more details see RB06). 



significantly in the viscous cases. Consequently, the stripped 
gas remains in larger clumps and is less readily mixed with 
the ambient medium. 

Also the width of the tail, or the flaring angle, is slightly 
smaller in the viscous runs. In the inviscid run, the maximum 
width of the tail is « 2 x 40 kpc, while it is ~ 2 x 35 kpc in all 
viscous runs. Thus, it is also unlikely that the small width 
of t he tail in NGC 4388 and the X-ray tails of ES0137- 
001 |Sun et alj|2006l ) can be explained by a microphysical, 
Spitzer-type viscosity. Interestingly, the tail appears to be 
shorter, the larger the viscosity. This is due to the fact that 
larger gas clouds, even if they have the same density like 
smaller clouds, experience a smaller acceleration for a fixed 
ram pressure. We note that the mass fraction of stripped gas 
in dense form does not depend strongly on viscosity. In the 
viscous cases, there are a few large dense clouds, whereas in 
the non- viscous case there are numerous small dense clouds. 
We do not show any plots here for the supersonic runs but 
they show the same qualitative behaviour. 

Figure [5] compares the evolution of the mass and ra- 
dius of the remaining gas disc for viscous and inviscid runs. 
The top panel demonstrates that the viscosity has nearly 
no influence on the evolution of the gas disc's radius. The 
bottom panel displays the gas mass in a cylinder of radius 
27 kpc and a height of 10 kpc. Interestingly, the mass evolu- 
tion in all runs is very similar, for all values of the viscosity 
that we have explored. The first dip in the graph marks the 
end of the instantaneous stripping phase. Very slight dif- 
ferences introduced by the viscosity occur during the next, 
the intermediate phase and also during the final continuous 
stripping phase. We have also measured the gas mass bound 
to the gravitational potential of the galaxy and its behaviour 
is shown in the middle panel of Fig. [2] The differences be- 
tween runs with and without viscosity is similar to that of 
the gas disc mass. This implies that the amount of gas lost 
from the galaxy is fixed by the ram pressure. 

This result shows that the viscosity in the parame- 
ter range considered here has a minor impact on the mass 
loss history of the gas disc and that ram pressure push- 
ing is the dominant mechanism, as f ound in earlier papers 
fe.g. lRoediger fc Bruggenll2006l . l2007l ). The fact that the vis- 
cosity has a minor impact on the mass loss history should 
also hold for most galaxies which move though clusters and 
thus experience a variable ram pressure - for the very rea- 
son that ram pressure pus hing is the main cause for mass 
loss also for these galaxies l)Roediger fc Briiggcn 2007). The 
mass loss due to ram pressure pushing, i.e. during the first, 
instantaneous stripping phase, is especially insensitive to dif- 
ferent ICM viscosities. An exception may occur when galax- 
ies move near edge- on. 



3 RESULTS 

Figure[T]shows the density in the subsonic runs with increas- 
ing viscosity: /„ = 0,0.1 and 0.5. The different behaviour of 
the galactic wakes is evident: The larger the viscosity, the 
less turbulent the wake. In the viscous cases, the ICM flows 
past the galaxy more smoothly than in the non-viscous runs. 
In the non-viscous cases, the KH- and RT-instability cause 
vortices and turbulence. Both instablilities are suppressed 



4 DISCUSSION 

Clearly, we have to be careful with quantitative predictions 
especially concerning the morphology of the wake because 
our simulations are two-dimensional. However, in previous 
work we have shown that mass loss rates are the same in 2D 
and 3D simulations. Moreover, in the non-viscous cases, the 
wake structures in 2D and 3D were very similar. Thus, our 
main conclusions should be robust: 
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Figure 1. Contour plot of the density in subsonic runs at 300 Myr, for increasing viscosity from left to right: left is inviscid (/ v = 0), 
middle is viscous with f v = 0.1, right is viscous with f v = 0.5. 



• The mass loss from the gas discs is hardly influenced 
by viscosity. 

• With increasing viscosity, the wake shows less stucture 
and turbulence, but larger clumps. 



A detailed investigation of the fate of the stripped 
gas requires several additions to our simplified model: a 
3D treatment to allow for other inclinations than face-on, 
prescriptions for heating, cooling and thermal conduction. 
These processes do not only influence the temperature in 
the stripped gas, but also determine in which wavebands the 
stripped gas will be observable. For the fate of the stripped 
gas, viscosity may make a difference as it may affect the 
thermal history of the stripped gas. Gradients of density and 
temperature are smeared out in the wake in the presence of 
viscosity. This will affect the efficiency of heat conduction 
and evaporation of clumps in the galactic wakes. It will also 
affect the efficiency with which mater ial in the wake can form 
stars as observed bv lSun et aLl(|2007n . Comparisons between 
models and observations will reveal which processes are the 
dominant ones in shaping galactic tails. 
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